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ABSTRACT 



Q0151+048 is a physical quasar (QSO) pair at z ~ 1.929 with a separation of 3.3 arcsec on the sky. In the spectrum of the brighter 
member of this pair, Q0151+048A, a damped Lya absorber (DLA) is observed at a higher redshift. We have previously detected the 
host galaxies of both QSOs as well as a Lya blob whose emission surrounding Q0151+048A extends over 5 x 3.3 arcsec. We seek to 
constrain the geometry of the system and understand the possible relations between the DLA, the Lya blob and the two QSOs. We also 
aim at characterizing the former two objects in more detail. To study the nature of the Lya blob, we have performed low-resolution 
long-slit spectroscopy with the slit aligned with the extended emission. We have also observed the whole system using the medium- 
resolution VLT/X-shooter spectrograph and the slit aligned with the two QSOs. The systemic redshift of both QSOs is determined 
from rest-frame optical emission lines redshifted into the NIR. We employ line-profile fitting techniques, to measure metallicities and 
the velocity width of low-ionization metal absorption lines associated to the DLA, and photo-ionization modeling, to characterize the 



DLA further. We measure systemic redshifts of : 



1.92924 + 0.00036 and 



<.mi(B) 



1.92863 ± 0.00042 from the H/3 and Ha 



emission lines, respectively. In other words, the two QSOs have identical redshifts within 2cr. From the width of Balmer emission 
lines and the strength of the rest-frame optical continuum, we estimate the masses of the black holes of the two QSOs to be 10'"'^ 
Mo and lO**"^'* Mq for Q 0151 +048A and Q 0151 -H048B, respectively. We then use the correlation between black hole mass and dark 
matter halo mass to infer the mass of the dark matter halos hosting the two QSOs; lO'^ ^"* Mq and 10" ''' Mq for Q 0151 -I-048A and 
Q0151+048B, respectively. We observe a velocity gradient along the major axis of the Lya blob consistent with the rotation curve of 
a large disk galaxy, but it may also be caused by gas inflow or outflow. We detect residual continuum in the DLA trough which we 
interpret as emission from the host galaxy of Q0151-I-048A. The derived H'' column density of the DLA is logW^o = 20.34 + 0.02 
cm"^. Metal column densities are also determined for a number of low-ionization species resulting in an overall metallicity of 0.01 Zq. 
We detect C ll* which allows us to make a physical model of the DLA cloud. From the systemic redshifts of the QSOs, we conclude 
that the Lya blob is associated with Q0151+048A rather than with the DLA. The DLA must be located in front of both the Lya blob 
and Q0151-I-048A at a distance larger than 30 kpc and has a velocity relative to the blob of 640 + 70 km/s. The two quasars accrete 
at normal Eddington ratios. The DM halo of this double quasar will grow to the mass of our local super-cluster at z = 0. We point out 
that those objects therefore form an ideal laboratory to study the physical interactions in a z = 2 pre-cursor of our local super-cluster. 
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1. Introduction 

Searches for Lyo' emission in the high redshift Universe have 
in the l ast decade found many sources of a type called Lya 
"blobs" jFvnbo et alJll999l:lKeel etal.ll 19991: ISteidel et al.l|200 " 



[Francis et al.' '2001; Matsuda et al.' '2004'; 'Palunas et al.' '200 
Weidinger et al.,.2004; Dev et al...2005; Villar-Martin et al..2005c 



Send offprint requests to: tayyaba@dark-cosmology.dk 

* Based on observations done with i) European Southern Observatory 
(ESO) utilizing 8.2m Very Large Telescope (VLT) X-shooter spec- 
trograph on Cerro Paranal in the Atacama Desert, northern Chile, ii) 
2.56m Nordic Optical Telescope (NOT), a scientific association be- 
tween Denmark, Finland, Iceland, Norway and Sweden, operated at 
Observatorio del Roque de Los Muchachos on the island of La Palma, 
Spain. 



iNilsson et al.ll2006h . Lya blobs must be considered the extreme 
end of the distribution of Lya emitting sources, in that they are 
very extended (sizes more than several times ten kpc across), and 
unusually bright in the Lya line (log Lya > 43 erg s"'). They are 
interesting and promising cosmological probes, partly because 
they are some of the largest cohesive structures at high redshift 
and partly because they hold key insights into the energetics of 
their source. Among the Lya blobs, a number of sub-categories 
have emerged, sorting them based on the presumed main source 
of ionizing flux. Several bl obs have been shown to be powered 



by vigorous star formati on ( Taniguchi et al.l200lt lOhvama et al.l 
120031 iMori et al.1 l2004t iMatsuda et alj l2006h). others by cold 
gas falling onto a da rk matter (DM) halo ( Nilsson et al J 120061 ; 
ISmith & Jarvisll2007h . 
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In this group of extreme objects, the most extreme sub- 
sample in terms of sizes and fluxes are those powered by ac- 
tive galactic nuclei (AGN), QSOs and radio galaxies. Gigantic 
Lya haloes a re routinely discovered aro und high redshift radio 
galaxies (see iMilev & De Breuc3 (l2008h for a revi ew). Smaller 
blobs are also often powered by AGN activity; iGeach et"al] 
(l2009l) find that roughly 17% of narrow-band selected Lya 
blobs are powered by AGN. In especially rare cases, multiple 
component systems including a QSO, a DLA, and a Lya blob 
are found (iM^ller et al.l [19981 iLeibundgut & RobertsonI Il999t 
iHennawi et al. I l2009h . It is yet unclear how these systems inter- 
act, or even their actual spatial configurations. 

One such system is the binary QSO Q0151-H048A & B 
(also called PHL1222), who se brightes t mem ber was first ob- 
served spectroscopically by iBurbidgj (1 19681) who measured 
its redshift to be z = 1.943 (at that ti me the other member 
of the binary was not yet discovered). IWilliams & WevmannI 
([1971 were the first to discuss the nature of the Zahs > 
Zem DLA in the spectrum of PHL1222 in some detail. They 
inferred a redshift of the absorber of more than 2000 km 
s"' relative to the redshift of Q0151+048A (hereafter qA). 



Table 1. Log of spectroscopic Observations of Q 015 1 +048 A & 
B with NOT (ALFOSC) and VLT/X-Shooter (XSH). 



The DLA is hence a member 
DLAs (see lM0ller et al 



199 



of the class of Zabs ~ Zem 

which later be came proximate 



DLAs (PDLAs) jEUison et all (l2010t) : see also IWevmann et alj 
(Il977h :l Prochaska et al.l (l2008bh ) Compared to normal interven- 
ing DLAs PDLAs tend to have evidence for a higher radiation 
field (presumably t he proximate QSO) and on average slightly 
higher metallicities dElhson et al 

mm- 

The QSO redshift reported in the above work was measured 
from the UV emission lines from N"^^, Si^^ and C"*^, which 
are now known to be systematically b lu eshifted relative to the 
system ic redshift (iTvtler & Fad Il992h . IWilliams & WevmannI 
d 19761) also inferred a minimum distance of 0.4 Mpc between 
the DLA and qA base d on the absence of fine structure lines 
from Si^ and C^. Later iMevlan et all (Il990l) were the first to es- 
tablish that PHL1222 has a fainter companion at ne arly the same 
redshi ft. The pair was discussed in more detail in lM0ller et alJ 
(119981) where Lyg emission was also detected in the DLA 



trough. iFvnbo et al . (1999, hereafter FMW99) found using nar- 
row band imaging that there is very extended Lyo' emission close 
to the redshift of the DLA and in retrospect this appears to be one 
of the first detections of what is now known as Lya blobs. 

In this study we present new spectroscopic observations of 
both members of the QSO pair Our main objectives are i) to 
characterize the DLA absorber (in particular its metallicity and 
velocity profile), ii) to determine the systemic redshift of both 
members of the pair using rest frame optical emission lines. Hi) 
further constrain the geometry of the system by looking for ab- 
sorption in the spectrum of Q 015 1 +048B (qB) at the redshift of 
the DLA, and iv) to further study the nature of the Lya emission 
blob. 

The outline of the paper is as follows: In §2 we describe our 
spectroscopic observations and data reduction. In §3 & 4 we an- 
alyze emission features (the Lya emitting gas and QSO emission 
lines) and absorption signatures of the DLA respectively and in 
§5 we discuss how these new data change our interpretation of 
the geometry and dynamical state of the system. In §6 we briefly 
summarize our conclusions. Throughout this p aper, we assume 



a Hub ble constant of Hq - 72 km s Mpc 



dFreedman et al.l 

I2OOI and a universe with flat cosmology Q,„ = 0.27, Qa = 0.73 
following the results presented in the initial WMAP data release 
(Spergel et al. 2003). At a redshift of 1.93 and with this cosmol- 
ogy, 1 arcsec corresponds to a proper length of 8.384 kpc. All 
wavelengths are corrected to vacuum heliocentric frame. 



Date 


Object 


Exp. 


Instrument 


slit 


R 


PA 






(sec) 




width 




(deg) 


2 Sep, 97 


A+S4 


14000 


ALFOSC 


0.7" 


700 


-78.5 


3 Sep, 97 


A+S4 


15000 


ALFOSC 


0.7" 


700 


-78.5 


4 Sep, 97 


A+S4 


4000 


ALFOSC 


0.7" 


700 


-78.5 


18 Nov, 08 


A+B 


3600 


XSHAJVB^ 


0.8" 


6200 


41.2 




A+B 


3600 


XSHA'ISJ 


0.7" 


11000 


41.2 


29 Sept, 09 


A+B 


960 


XSHAJVB 


1.0" 


5100 


41.2 




A+B 


960 


XSHA'IS 


0.9" 


8800 


41.2 




A+B 


480 


XSH/NIR 


0.9" 


5100 


41.2 



" UVB-arm range: 3000-5500 A 
* VIS-arm range: 5500-10000 A 



2. Observations and data reduction 

2.1. The NOT/ALFOSC data 

On three nights in September 1997 we used the NOT 
equipped with Andalucia Faint Object Spectrograph and Camera 
(ALFOSC) to obtain a spectrum of the elongated Lya emission 
line object discovered by FMW99 (see their Fig. 2) who named 
it S4. Our slit was centered on qA and aligned with the major 
axis of S4 at a position angle of -78.5 deg (see Fig.[T]i. We used 
ALFOSC grism G6 and a slit of 0.7" which provided a resolu- 
tion of R - A/ AA - 700. We integrated for a total exposure time 
of 33000 sec at this slit position. A detailed log is provided in 
Table [Tj The data were reduced using standard techniques for 
bias and flat-fielding and a final stack was obtained using opti- 
mal weights. 

2.2. VLT/X-shooter data 

On November 18th 2008, during the first instrument commis- 
sioning run, we returned to this field, this time using the X- 
shooter on the VLT which allowed us to observe both QSOs 
with higher spectral resolution and long wavelength cover- 
age from ultraviolet (UV) to near-infrared (NIR). The ESO- 
VLT X-shooter is an Echelle spectrograp h mounted at the VLT 
Cassegrain focus dD'Odorico et al. I I2OO6I) . Here we used a slit 
PA of 41.2 deg which placed both QSOs on the slit (see Fig.[T]) 
thereby allowing us to obtain a 3600 seconds exposure of both 
QSOs (Program ID: 60.A-9022(C)). During this early commis- 
sioning run only two of the three arms were operational and we 
obtained spectra in parallel with the UV-blue (UVB) and visual 
(VIS) arms of the spectrograph using slits of 0.8" and 0.7" re- 
spectively. The resulting resolutions are 6200 (UVB) and 1 1000 
(VIS). Because the seeing was significantly larger than the slit 
widths the sky-line widths provide a good measure of the reso- 
lution. 

On September 29, 2009, during a Science Verification run 
(Program ID: 60.A-9441(A)), we again integrated with the X- 
shooter slits on both QSOs. This time we could utilize all three 
arms (i.e. UVB, VIS, and NIR) and the log of observations is 
again provided in Table [T] Availability of the NIR arm during 
the September run allowed us to extend the wavelength coverage 
through to rest wavelength optical lines. We obtained 480 sec in 
the NIR arm with a resolution of R = 5100. All X-shooter data 
are available from the ESO archive. 

The X-shooter data were processed using a pr ehminary ver- 
sion of the X-shooter data reduction pipeline dGoldoni et al.l 
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Fig. 1. Top frame is a cut-out of the original narrow band image 
(FMW99). Contours are shown in black and the centres of the 
two QSOs are marked by a white cross. NE is to the upper left. 
In the bottom frame the QSOs were psf-subtracted and only the 
extended Lya emission is seen. The crosses still mark the po- 
sition of the QSOs. The slit-position along the Lya major axis 
used for the NOT/ALFOSC observations is shown in green, the 
X-shooter A-i-B slit-position is shown in red. 



l2006HModigliani et al ] l2010h . The pipeline first corrected UVB 
and VIS raw frames for bias and NIR frames for dark current. 
After division by the flat field, the pipeline performed back- 
ground subtra ction, corrected cosm ic ray hits using the method 
developed by van Dokkum (2001') and subtracted sky emission 
lines using the Kelson (2003) method. The individual orders 
were extracted and rectified in wavelength space. This was done 
by using the previously obtained wavelength solution from cal- 
ibration frames. The individual orders were merged afterwards 
and in the overlapping regions the merging was weighted by the 
errors which were being propagated during the reduction pro- 
cess. From the final 2D merged spectrum one dimensional spec- 
tra for qA & qB were extracted. 

We obtained a spectrum of the hot B-type main sequence star 
Hip02097 1 with the NIR arm on the same night. The spectrum 
was extracted using the same procedure and was subsequently 
used for atmospheric absorption correction. No flux calibration 
was applied to either of the spectra. 

3. Results 

3.1. The NOT/ALFOSC data 

The purpose of the NOT spectrum was to study the dynamics of 
the Lya emitting gas (S4). The resolution of the spectrum is too 
low to provide any useful information about the DLA metallicity, 
so the primary interest is in the region around Lya emission. 
In Fig. |2] (left panel) we show this region of the 2D stacked, 
sky subtracted spectrum. We here use a negative representation 
where black is bright and orange is zero. In the centre the DLA 
absorption trough is a dominant feature, and the arrow points to a 




Fig. 2. Left panel: NOT/ALFOSC 2D spectrum of qA, the dis- 
persion is vertical with short wavelengths at the top. The DLA 
line is at the centre, the arrow points to faint extended Lya emis- 
sion. Right panel: Same as left after SPSF subtraction. The full 
extent of the Lya emission is now visible. It is seen to be tilted 
and shifted relative to the DLA. The vertical axis provides rela- 
tive velocity (km/s) in the rest frame of the DLA absorber 
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Fig. 3. Velocity field of the Lya emission seen in the right panel 
of Fig. 12] The DLA absorption sightline towards the qA is 
marked by a " ★ ". 



small but significant area of extended line emission. To examine 
the Lya emission closer to the QSO it is necessary first to remove 
the emission from the QSO , a technique referred to as Spec tral- 
PSF (or SPSF) subtraction (lM0lleiil2()OOl: lM0ller et al.llIOOOh . 

In the right panel of Fig. |2] we show the same part of the 
spectrum after SPSF subtraction around the extended Lya emis- 
sion. We clearly detect extended Lya over a total region of about 
5 arcsec (42 kpc proper), corresponding well to the extend of the 
emission seen in narrowband imaging (Fig. [1] lower panel). It 
is also clear that the emission is not symmetrically distributed 
around the QSO (also consistent with Fig.[T]i, that it is well sep- 
arated (in redshift) from the DLA and, somewhat surprisingly, 
that the emission has a strong velocity gradient along the major 
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Fig. 4. X-shooter spectra of both QSOs obtained on November, 
18 2008 with a single exposure of 3600 sec. Upper panel: 2D 
UVB-arm spectra of both QSOs, qB is the fainter object and qA 
is the brighter object displaying strong DLA absorption. Middle 
panel: ID spectrum of qB. Lower panel: ID spectrum of qA 
showing the DLA absorption. Note that in both the 2D and the 
ID spectrum of qA residual flux is clearly visible in the damped 
trough. 



axis which has a remarkable resemblance to a rotation curve. The 
velocity gradient as measured on the SPSF subtracted spectrum 
is shown in Fig.|3] 

3.2. VLT/X-shooter data 

3.2.1 . Residual flux in the DLA trough 

One of the most striking features of the X-shooter spec- 
trum of qA is that the DLA absorption line does not 
go completely to zero. Lyo- emission in or close to 
DLA absorption lines has been s een in several cases. 
lM0ller & Warred ([19931). iPettini et alj (Il995h. iDiorgovski et al.l 



(119961). iLeibundgut & RobertsonI (Il999h. lM0ller et al.l (l2002h . 
lM0ller et al.l ( l2004l) and lFvnbo et al.l ( l2010l) all reported a narrow 
emission line slightly off'set from the position of the background 
quasar and consistent with being Lyo' emission fr om the DLA 
host w hich absorb the QSO to zero flux at Lyo-. ICooke et al.l 
(1201 Ol) report what they interpret as either a broad absorption 
line with emission at the centre or as two absorption lines sepa- 
rated by non-absorbed quasar continuum. Unfortunately they do 
not provide the detailed spatial information which would help 
the interpretation, but since a possible interpretation of the "fea- 
ture" is that it is non-absorbed quasar continuum it must be cen- 
tered on the position of the quasar and have the spatial profile 
of a point-source. This makes it unlikely that it is Lyor emission 
as Lyo- emission usually is more extended than a point-source 
and in addition it would be an unlikely coincidence that the Lyo- 
emission by chance would be lined up exactly with the second 
im age of the lens where th e DLA covers only the first image. 

iHennawi et al.l (l2009h also report an emission line in a 
damped absorption line. They do not provide any information 
about spatial offset, b ut the object has a large (5 arcsec) size 
along the slit. Henna wi et al ] (l2009h rule out that the emission is 
related to the DLA and instead interpret it as a Lya blob phys- 
ically connected to the quasar This object has two features in 



common with the DLA in qA. Both are PDLAs and both have 
Lyo' emission extending 5 arcsec. The extended Lyo- of qA is 
however significantly offset in redshift from the PDLA and does 
not enter into the absorption line. Nevertheless in both the 2D 
and the ID X-shooter spectrum of qA (Fig. |4]i it is clearly seen 
that there is a smooth continuum providing a non-zero "floor" in 
the Lyo- line. In this section we aim to determine the source of 
this residual flux in the damped absorption line. 

There are two possible ways to explain this residual flux. 
Either the absorber is so small that it only partly covers the emit- 
ter (the central engine), or there are at least two different sources 
of emission and the DLA absorber is only covering one of those. 
In Fvnbo et al.. t2000i) it was reported that qA is located inside 
a comparatively bright host galaxy with a B(AB) magnitude of 
20.8+0.2 and /(AB) of 21.5+0.2. The host was best fitted by a de 
Vaucouleurs profile, and in Z?(AB) it was found to be 3.0 magni- 
tude fainter than qA while in m(AB) and /(AB) it was 3.9 and 4.0 
magnitude fainter respectively. We now investigate if the resid- 
ual flux in the DLA trough simply is the spectrum of the QSO 
host galaxy. 

We first co-added the profile of all wavelength bins in the 
centre of the DLA line where it should be completely absorbed. 
We then did the same where the QSO continuum had recovered 
both on the blue side (3540-3545 A) and on the red side (3595 - 
3605 A). The FWHM of the QSO profile along the slit on both 
the blue and the red side is 1.03 arcsec while the FWHM of the 
flux at the bottom of the DLA is 1.57 arcsec, i.e. significantly 
wider. Simple subtraction in squares gives an intrinsic FWHM of 
1.2 arcsec for the object seen in the DLA line which immediately 
rules out that the residual could be caused by partial coverage of 
the QSO alone because in this case the profile should be the same 
as that of the QSO. The object that we see in the DLA trough is 
therefore a different object than the QSO, it is relatively large 
(1.2 arcsec FWHM) and could possibly be the host galaxy. 

The spectrum was taken with a slit of 0.8" and we therefore 
have slit losses of both the QSO and of the other object. We have 
modeled the psf of the QSO such that the psf-on-slit is identi- 
cal to that observed. Using this we measure a slit loss of 40%. 
Similarly we have modeled the profile of the object seen in the 
trough, and find slit losses of 60% in this case. The measured 
counts per bin of the extended object is a factor 27 less than 
of the QSO continuum. Correcting for the different slit losses 
this becomes a factor 18 fainter, i.e. 3.1 magnitudes fainter This 
corresponds well to the reporte d magnitude difference as found 
via imaging dFvnbo et al.ll2000() . From Fig.|4]upper panel (more 
clearly seen in Fig. |7ll it is also seen that the two objects are 
centered on the same position. 

Summing up all of the evidence, extended profile, centered 
on the QSO, similar brightness relative to qA, we conclude that 
the flux seen in the DLA trough is the signa ture of the object 
reported to be the host in dFvnbo et alj|2000h . It remains to be 
determined whether the flux from the host is stellar emission or 
dust reflection. The gradient of the emission in the trough could 
suggest that it at least partly is quasar Ught reflected off dust 
in the host halo. Spectro-polarimetric data would be required to 
resolve this question. 

3.2.2. Emission line analysis of Q0151+048A & B 

The original interpretation of the Lya emission from this sys- 
tem was severely complicated by the large apparent blueshift of 
qA relative to the DLA (1250 km s"') as well as to t he Lya 
emitt e r. Such shifts are common in QSOs spectra (see i Grandil 
119821: IWilkesI 119861: lM0ller & Warred fl993E iLaor et alJ 1199.51: 
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Table 2. Emission line properties of Q0151+048A & B. The columns provide observed and restframe wavelengths, observed and 
restframe equivalent widths with their corresponding errors, line IDs and emission redshifts. 



Object 








C" EW 




ID 






A 


A 


A 


A 


A 








qA: 


3612.1710 


1240.13 


13.658 


0.237 


4.688 


Nv 


1.91274 


± 0.00014 




3694.1012 


1262.59 


3.9000 


0.090 


1.330 


Sin 


1.92581 


± 0.00027 




3814.8126 


1304.35 


11.400 


0.113 


3.900 


Oi 


1.92468 


± 0.00035 




3894.1912 


1334.50 


3.8000 


0.087 


1.300 


Cll 


1.91809 


± 0.00016 




4063.6774 


1399.70 


40.350 


0.198 


13.90 


Si IV + Oiv 


1.90325 


± 0.00030 




4504.6833 


1549.06 


62.700 


0.266 


21.56 


Civ 


1.90801 


± 0.00040 




6804.3376 


UV 344(2325.38) 


19.160 


0.110 


6.550 


Fen 


1.92612 


± 0.00026 




7109.3696 


UV 360(2433.436) 


10.730 


0.140 


3.670 


Fen 


1.92154 


± 0.00038 




8188.7607 


2798.75 


96.440 


0.240 


32.96 


Mgll 


1.92586 


± 0.00040 




14244.000 


4862.68 


119.80 


0.280 


40.90 


HyS 


1.92924 


± 0.00036 




19231.400 


6564.61 


525.56 


0.810 


179.4 


Hey 


1.92939 


± 0.00043 


qB: 


3557.3960 


1215.67 


243.70 


0.923 


83.300 


Lya 


1.92628 


± 0.00034 




3628.3042 


1240.14 


8.2900 


0.631 


2.8300 


Nv 


1.92572 


± 0.00042 




3678.1892 


1262.59 


2.1300 


0.409 


0.7300 


Sin 


1.91321 


± 0.00045 




3810.1611 


1304.35 


2.3840 


1.009 


0.8190 


Oi 


1.92112 


± 0.00036 




4094.7356 


1399.70 


32.023 


0.963 


10.986 


Si IV + Oiv 


1.92544 


± 0.00025 




4528.4208 


1549.06 


60.699 


1.039 


20.752 


Civ 


1.92333 


± 0.00043 




8198.4067 


2798.75 


89.800 


0.270 


30.650 


Mgii 


1.92931 


± 0.00039 




19225.300 


6564.61 


492.72 


2.580 


168.20 


Hq- 


1.92863 


± 0.00042 





Q0151+048A 
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Fig. 5. Ha emission lines from Q 015 1 +048 A & B after subtrac- 
tion of the QSOs continuum. As described in the text only the 
good data intervals between the strong airglow residuals were 
used. The data, including airglow residuals are shown as a light 
dotted line, the solid black line is data with the not used parts 
removed. The solid red Unes show the best gaussian profile fit 
to the used intervals. The inferred redshifts are ZHa{A) = 1.9294 
and ZHa(B) = 1.9286. 



Mcintosh etaLl 119991: IVanden Berket al.l 120011: iRichards et all 
2002h . but even though there is a general pattern to the shifts 
it has not yet been possible to device a method to successfully 
correct for this. We therefore had serious doubts about the accu- 
racy of QSO systemic redshifts determined from UV emission 
lines, and one of the two main drivers for the X-shooter spec- 
troscopy was to obtain systemic redshifts based on rest-frame 
optical emission lines. This is possible because of the X-shooter 
NIR arm. In Fig.|5]and|6]we show the H a and H/? emission lines 
respectively. For qB, the fainter of the two QSOs, only H a was 
detected. The blue part of the H a emission lines was affected by 
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Fig. 6. Hj6 emission line from Q 015 1+048 A after the contin- 
uum subtraction. The solid red line shows the best gaussian fit 
(see caption to Fig.|5]for details) providing an emission redshift 
of ZHAA) = 1.9292. 



atmospheric absorption which corrected out well with the use of 
hot B-type main sequence stars. 

In Table |2] we list observed and rest-frame parameters for 
the emission lines of the spectra of qA & qB. The continuum is 
well detected for both QSOs and the equivalent widths are mea- 
sured by integrating the emission line between the end points and 
above the estimated continuum. For the Anh s and inferred z^v w e 
use the same peak (or mode) definition as in| M0lleret al.l([l998l) . 
It is seen that the inferred redshifts of qA range from 1 .9034 to 
1 .92939 corresponding to a relative velocity of 2700 km s"' vin- 
dicating our doubts about the previous redshift determination. 

The narrow and unblended [O ill] A5QQ7 emission line is pre- 
sumed to trace the systemic (center of mass) red shift of the 
galax y ^Gaskell 1982; Vrtilek & Carleton 1985; M cintosh et al.l 
Il999h . Unfortunately the [O ill] emission lines at our redshift are 
in a region of strong telluric absorption. The telluric absorption 
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combined with the very narrow profile of [O ill] makes it impos- 
sible to recover the lines via telluric correction. The H a line is 
in a part of the spectrum which is even stronger absorbed, but the 
width of the line is here helpful in recovering parts of the line. 
The H j6 line is intrinsically fainter than the H a but in a region 
of less telluric absorption and therefore easier to recover. 

In addition to the telluric absorption there are many strong 
sky emission lines present causing large residuals after sky sub- 
traction. With the high resolution of X-shooter and for the wide 
Ha and H/3 lines this is not problematic as we simply masked 
out the affected intervals and used only the remaining sections 
(see Fig. |5]l. We subtracted continuum from the emission lines 
and then proceeded to fit a smooth profile to the central part of 
the lines using only the intervals not masked out. After the fit to 
the data we used the Icr error spectrum to create 1,000 random 
realizations of errors, added those to the spectrum and repeated 
the process. In Table |2] we list the standard deviation of this 
distribution as statistical errors on the redshift. For qA we com- 
pute the optimally combined redshift as the weighted mean of 
the two Unes and obtain Zem(A) - 1.92930 ± 0.00028, for qB we 
only have Ha and therefore simply adopt the redshift of that 
line as the systemic Zem(B) - 1.92863 + 0.00042. Redshifts de- 
termined from those lines are known to have an offset of only a 
few tens of km s ^ from the true systemic redshift of the QSO 
(IVanden Berk et al.ll200lh . 

3.2.3. Extended emission in tine X-sliooter spectrum 

In Fig.[T]it is seen that the X-shooter slit follows a line across the 
Lya emitter more closely aligned with it's minor axis. We would 
still expect to see some extended emission, but not extending as 
far away as for the NOT spectrum. To confirm this expectation 
we therefore also performed SPSF subtraction on the X-shooter 
spectrum. The result is shown in Fig.|7] where we have marked 
the redshifts of both QSOs relative to that of the DLA. SPSF sub- 
traction of qA was only performed in the range -1300 to -200 
km s ' (relative to the DLA redshift). This procedure left the flux 
from the host galaxy inside the DLA trough and clearly shows 
that the damped line does not at any point go to zero flux. Also 
it is seen that the host is precisely lined up with the position of 
qA. The extended Lya is also here seen to be asymmetric and 
to stretch about half way towards qB, consistent with the NB 
imaging (Fig.[T]i. The emission is here offset 750 km s"' from 
the DLA, it has a FWHM of 200 km s ' and there is no obvi- 
ous evidence for a rotation curve along this axis. The difference 
in velocity offset at the position of the quasar seen at the two 
different slit PAs is due to the integration of the velocity field 
over different parts of the blob. The velocity measured along the 
major axis is less affected by this and we obtain a weighted com- 
bined offset of 640 + 70 km s"' . 

There seems to be marginal evidence for a double-emission 
profile with a splitting of 140 km s"' but the S/N of the 
data is too low to make any conclusive statements about this. 
D ouble-pe aked emi ss ion line profiles from Lya were predicted 
bv lHarrin gton'('l973') from a static, un-absorbing, opaque scat- 
tering medium. Verhamme et al. (2006) obtained the same re- 
sult using a 3D Lya radiative transfer code with a Monte- 
Carlo fitting method, and subsequently used the code to fit Lya 
lines in a sample of hi gh redshift Lyman-break galaxies (LBGs; 
IVerhamme et al.ll2008h . In this publication, two LBGs had dou- 
ble emission line peaks and were best fitted with a dust-free, 
nearly static medium. Thus, if it can be confirmed that the emis- 
sion line profile is indeed double-peaked then this could be indi- 



z(em,B) 

I 




Fig. 7. 2D X-shooter spectrum with the slit placed on both 
QSOs. Short wavelengths are to the left. The centre of the DLA 
line is marked and chosen as fiducial zero velocity. SPSF sub- 
traction was done for the qA spectrum in the range -1300 to 
-200 km s"'. There is a clear extended residual seen which we 
interpret as Lya emission at -750 km s"'. The redshifts of both 
QSOs based on H a and H/? are also marked. 



eating that the blob is static rather than infalling or the signature 
of an expanding wind. 



4. Absorption line analysis 

The X-shooter absorption line spectra are plotted in Fig.[TT] The 
medium resolution of X-shooter allows us to perform absorp- 
tion line profile fitting and determine column densities from lines 
that are not heavily saturated. Our second main driver for the X- 
shooter spectroscopy was therefore to derive ionic column den- 
sities and element abundances for the DLA system observed in 
the spectrum of qA. 

Element abundances are important for our understanding of 
this system, first because it has been suggested that the DLA 
metallici ty is related to the mass and lurn i nosity of the DLA 
galaxy (M0ller et al.' '2004^; iLedoux et al] l200a iFvnbo et all 
2DQ8; Prochas ka et_al, 2i)08b|) second because it can provide us 
with limits on the QSO/DLA distance 04.2l i which are needed in 
order to determine the arrangement of the QSO/DLA/Lya emit- 
ter system in 3D space rather than redshift space. In addition 
several systems at lower redshifts are also seen in the X-shooter 
spectra and for completeness we also identify those ( 94.31 ). 

4. 1 . Absorption lines of tiie DLA system 

lM0ller et all (1 19981) reported an H" column density of the DLA 
of log A^H" = 20.36. In our higher quality data we determine a 
column of log A^H" = 20.34+0.02 cm^^ at z = 1.93421+0.00001. 
The Voigt profile fit to the DLA is overplotted in blue on the 
normalized spectrum in Fig. [TT] At the same redshift we also 
identify several metal lines in the spectrum of qA, of which none 
are seen in the spectrum of qB down to a 3(t upper limit of 0.08 A 
(see Fig.[n]i. We identify lines from O", C+, C+* C^+, Si+, Si'^, 
Si^+, N°, Fe+, AL, Al^^, Mg° and Mg+. Equivalent widths and 
corresponding 1 sigma errors of those lines are listed in Table [3] 
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Table 3. Voigt-profile fitting (best fits shown in Fig. [8]) of metal ion transitions in the DLA using FITLYMAN. Ionic column den- 
sities were u sed to derive individual element abundances in the neutral gas phase via comparison to the solar neighbourhood 
dAsp lund et a l.l2009l) . Multiplets were fitted with same column densities and turbulent broadening parameter values. The Si ii /II 808, 
Fe II /il611 and N v lines are not detected and for those we provide 3cr column density upper Umits (see text). 



^obs 




O" EW 




ID 




logyv 


blurb 


[X/H] 


A 


A 


A 


A 








km s"' 




3634.960 








Nv(1238) 


1.93421 


< 12.71 






3646.650 








Nv(1242) 


1.93421 


< 12.71 






3519.770 


0.070 


0.030 




Ni(1199) 


1.93424 


12.95 ±0.12 






3820.900 


0.681 


0.029 


0.232 


1 (1302) 


1.93426 


> 14.97 


20.1 ±2.1 


[0/H] > -2.06 


3915.800 


0.724 


0.037 


0.247 


C II (1334) 


1.93421 


> 14.46 


18.8±0.2 


[C/H] > -2.31 


3919.150 


0.050 


0.015 


0.017 


Cii* (1335.7) 


1.93414 


12.94 ±0.12 






4542.600 


0.928 


0.040 


0.316 


Civ (1548) 


1.93413 


> 14.56 


18.9±0.4 




4550.155 


0.813 


0.035 


0.277 


Civ (1550) 


1.93413 


> 14.56 


18.9±0.4 




3492.924 


0.790 


0.023 




Sill (1190) 


1.93421 


a 






3501.357 


0.362 


0.014 


0.123 


Sin (1193) 


1.93421 


13.92 ±"0.03 


18.6±0.4 




3698.343 


0.767 


0.034 


0.261 


Sill (1260) 


1.93421 


13.92 ±0.03 


18.6±0.4 




3827.300 


0.265 


0.030 


0.100 


Sill (1304) 


1.93421 


13.92 ±0.03 


18.6±0.4 




4479.678 


0.459 


0.026 


0.156 


Sin (1526) 


1.93421 


13.92 ±0.03 


18.6±0.4 


[Si/H] = -1.93 ±0.04 


5305.090 








Sill (1808) 


1.93421 


< 14.46 






3540.130 


0.804 


0.030 


0.270 


Sim (1206) 


1.93421 


> 13.85 


18.8±0.3 




4089.480 


0.640 


0.065 


0.218 


Si IV (1393) 


1.93415 


> 13.75 


19.2±0.7 




4115.933 


0.517 


0.035 


0.180 


Si IV (1402) 


1.93415 


> 13.75 


19.2±0.7 




3359.492 


0.103 


0.024 


0.035 


Fell (1144) 


1.93421 


13.60 ±0.02 


18.8±0.3 


[Fe/H] = -2.24 ±0.03 


4719.540 


0.144 


0.035 


0.049 


Fell (1608) 


1.93421 


13.60 ±0.02 


18.8±0.3 




4727.600 








Fell (1611) 


1.93421 


< 14.67 






6967.170 


0.138 


0.010 


0.047 


Fe II (2374) 


1.93421 


13.60 ±0.02 


18.8±0.3 




6991.530 


0.883 


0.010 


0.301 


Fell (2382) 


1.93421 


13.60 ±0.02 


18.8±0.3 




7589.770 


0.410 


0.013 


0.140 


Fell (2586) 


1.93421 


13.60 ±0.02 


18.8±0.3 




4902.464 


0.321 


0.041 


0.110 


Alll (1670) 


1.93423 


12.46 ± 0.02 


18.8±0.2 


[Al/H] = -2.33 ± 0.03 


5442.160 


0.09 


0.023 


0.030 


Aim (1854) 


1.93423 


12.15 ±0.07 


18.8 ±0.2 




5465.850 


0.060 


0.030 


0.020 


Aim (1862) 


1.93423 


12.15 ±0.07 


18.8 ±0.2 




8205.080 


1.611 


0.010 


0.550 


Mg II (2796) 


1.93421 


> 13.79 


18.8±0.3 


[Mg/H] > -2.15 


8226.140 


1.380 


0.011 


0.470 


Mgii (2803) 


1.93421 


> 13.79 


18.8±0.3 




8371.210 


0.058 


0.021 


0.020 


Mgi(2852) 


1.93421 


11.19±0.11 







" The Si II /II 190 line is blended with a Lya forest line and was not 
included in the overall fit. 



We u sed the Voigt profile fitti ng FITLYMAN package in 
MIDAS dFontana & Ballestej [19951) and adopted the list of 
atom ic data, labora tory wavelengths and oscillator strengths 
from lMortonI (l2003h . FITLYMAN finds the best global fit using 
a x~ minimization algorithm which includes the spectral reso- 
lution and returns best fit parameters for redshift, column den- 
sity and Doppler turbulent broadening as well as errors on each 
quantity. It may be used to fit single lines one at a time, or one 
may (optionally) choose to fit several lines with the same param- 
eters (one may e.g. insist that all the chosen lines have the same 
redshift or turbulent broadening). 

We employed FITLYMAN in two ways. First we chose to fit 
all lines with the same redshift and turbulent broadening. The 
resulting fits were mostly very good, but for a few atoms and 
ions (e.g. 0° and Si"*^) the fits were slightly offset from the data. 
We therefore chose to fit each ionic species independently al- 
lowing them to have different redshifts and turbulent broaden- 
ing. The resulting column densities did not change significantly 
from one procedure to the other so our results are not depend- 
ing on this choice. The FITLYMAN output parameters are listed 
in Table [3] and the resulting fits (red) are shown overlaid on the 
data in Fig. [8] Note that the Si II /11190 line is in the Lyman forest 
and clearly blended with an intervening line. That line was not 
used for the fit, but the predicted line is plotted for completeness. 



While the signal-to-noise ratio of our X-shooter spectrum 
is excellent, its resolving power is not high enough to resolve 
the lines which at this resolution may have hidden velocity sub- 
structure. Accurate column densities can only be derived from 
weak lines while the column densities of stronger lines may 
be underestimated due to hidden saturation. For this reason, we 
consider our fit results for O I, C II and Mg II to be lower limits. 
For Si II, Al II, Fe ii, C Ii*, Mg i and N i, hidden saturation could 
be an issue in case additional components with b values narrower 
than 3 km s"' are present. The non-detection of Si 11/11808 sets 
an upper Umit of lO'" * on the column den sity of any component 
with a of 3 km s"^ . lEllison et alj ( 1201 Oh presented a spectrum 
of qA of high resolution (^ 6 km s"' FWHM) with similar de- 
tection limits as the X-shooter spectrum. By way of comparison, 
we note that our column densities of Si^, Al^, Fe^, N" and C^* 
are 0.09 , 0.11, 0.10, 0.11 an d 0.06 dex lower than those mea- 
sured bv lEUison etal] ( l2010l) which we take as the error induced 
by the uncertainty on our derived b value. This shows that nar- 
row hidden components account for less than 0.1 dex of the total 
column densities of these ions. 



For non-detected lines, we report in Table[3]3cr column den- 
sity upper limits. For those, we computed the column densities 
assuming them to be on the linear part of the curve-of-growth 
rather than from the Voigt-profile fitting. 
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Fig. 8. Voigt-profile fits (red overlay) to the metal absorption lines of the DLA in the spectrum of qA. The derived parameters are 
provided in Table|3] The vertical dashed black line is our adopted zero velocity corresponding to z = 1.93421. For Siii/11190, we 
did not fit the line but show the expected line in blue. For N i, we only list the strongest line in the table but the fit was performed 
simultaneous to the triplet as shown. 



In the la st column of Table [3] we adopt the solar abun- 
dances from lAsplund et all (120091 to derive metallicities rela- 
tive to solar. We find a silicon abundance, which is commonly 
used to define DLA metalHcities, of [Si/ H] = -1.93 ± 0.04 
suggesting a metal-poo r system (see Fig. l9l lEedoux et aDl2006t 
iProchaska et al.ll2008al) . We determine a line velocity width of 
64 km s ' as defined in iLedoux et al.l (l2006h using Fe II/12382 
which, after correction for the instrumental resolution of 27 km 
s leads to 58 km s ' . In Fig.|9] we plot this toget her with the 35 
DLAs in the redshift range 1.7 to 2.43 taken from lLedoux et al.l 
(I2006h . 

One can see that the point follows the general metallic- 
ity/velocity width relation at z = 2, but that it falls within the 
10-15% lowest metallicities, and that it lies on the high AV side 
of the bisector fit (in the widest 15%). 

Collisions with electrons are the dominant excitation mech- 
anism to populate the upper fine-structure energy level of the 

ground state (see also Williams & Wevmann 197(|). The 
Cli*/}1335.7 line can therefore be used to set an upper limit 
on the particle density «ho- Unfortunately, the Cii line only 
provides a column density lower limit due to saturation, but 
assuming solar abundance ratio of Si/C we can use Si^ as a 
proxy and compute the column density. We find that : log 
Nc^./Nc* = -1.9, and from fig. 4 of ( ISilva & Viegasll2002[) . we 
determine a particle density of hh" = 25 cm"-'. 



4.2. Photo-ionization model of the DLA 

Based on the column densities reported above, we now make 
a simple photo-ionization model of the absorbing cloud. The 
main purpose of the model is to constrain the distance between 
qA and the DLA. For the model, we use Cloudy version 07.02 
described bv iFerland et al.l (Il998h . The Cloudy calculations are 
based on the following assumptions. The cloud is assumed to be 
a plane-parallel slab of gas with total neutral hydrogen column 
density log Nno = 20.34, metallicity 1/lOOth of solar, solar rela- 
tive abundance ratios, and a particle density in the neutral phase 
(«n») of 25 cm""*. The radiation field is composed of the Cosmic 
Microwave Background at z = 1.934 and, towards one side of 
the cloud, the radiation spectrum of a typical AGN with spe- 
cific luminosity at the Lyman limit Ly^ = 30.9 erg s"' Hz"'. The 
latter value has been determined using the flux-calibrated spec- 
trum obtained from NOT in August 2008 and a best-matched 
spectral slope of - 1.76 between rest-frame 850 and 1600 A (see 
iTelfer et al.ll2002l) . 

Using this model, we performed a series of photo-ionization 
equilibrium calculations for a range of QSO-DLA distances and 
computed ionic column densities. Keeping in mind that this is a 
simple model we can nevertheless apply our observational con- 
straints to get an estimate of the distance. 
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Fig. 9. Average metallicity vs. velocity width o f 35 DLA low 
ionization line profiles from lLedoux et al.) (l2006h . The solid line 
represents the linear least square bisector fit. The DLA from the 
present w ork is marked by the large diamond. For details on this 
figure see lLedoux et al.l (l2006l) . 

From log A^^i-* /^ai+ - 0.3 1 + 0.09, we find agreement for 
30+11 kpc while the upper limit log Nj^-i+/N-i^o < 0.14 sets a 
hard lower limit on the distance of 28 kpc. This model is there- 
fore in best agreement with a distance of 30-40 kpc assuming 
no local source of ionization. The column densities of the higher 
ionization gas (Si~^, Si^^ and C^^) are all lower limits but we 
note in passing that those limits are all consistent with the above 
distance estimate. 

4.3. Intervening absorbers 

The primary purpose of the X-shooter spectroscopy was to char- 
acterize the DLA and the QSO. In addition to the DLA four other 
C IV absorption systems at redshifts 1.2939, 1.3360, 1.5613, and 
1.6192 were identified in the spectra. All four C iv doublets are 
seen in both spectra (as seen in Fig. [TTI) indicating that C iv ab- 
sorbers at those redshifts typically have sizes larger than 30 kpc 
and covering factors close to unity. This is consiste nt with pre- 
vious studies (iRauch et al.ll200ll:lEllison et al.ll2004 . One of the 
absorbers (z - 1.6192) is also seen as Si iv absorber in the sight- 
line towards qA, but not towards qB, which might indicate that 
Si IV absorbers extend to smaller distances than C iv absorbers. 

It is striking in Fig.[TT]to see the one-to-one correspondence 
between the Unes in the Lyman forest of the two QSOs (see 
ISmette et al.lll992l Il995h . The few lines which are not repro- 
duced in the qB spectrum all belong to the DLA system. 

5. What is the Q 0151 +048A&B/DLA/Lya-blob 
system? 

This complex system has been the target of several studies be- 
cause it brings together several types of objects which at this 
redshift are usually seen only separately. The mix of objects has 
made the interpretation difficult, yet at the same time the interac- 



Table 4. BH, total masses and edingtton ratios of the QSOs. 
Masses are provided in logarithmic scale. 







qA 


qB 


r^ll 7T TTV if T T in /I — 1 \ 

FWHM HyS (km s ) 




3577 ± 150 




FWHM Mgll (km s"') 




2790 ±110 


1840 ± 130 


r WrIM IV (km S ) 




JO /j + 1 


o/i on _i_ 1 1 n 


ALa (10-*'' erg s"') 


(5 100 A) 


170 


5.8 


ALa (lO"* erg s"') 


(3000 A) 


2100 


105 


ALAiW** ergs-') 


(1450 A) 


1400 


67 


Mbh /Mq 


(H/J) 


9.12 




Mbh /Me 


(Mgll) 


9.37 


8.38 


Mbh I Me 


(Civ) 


9.50 


8.36 


M,JMe 




13.74 


13.13 


log i^bol/^-Edd 




-0.24 


-0.54 



tion between the objects allow us to extract information which is 
not normally accessible and can only be obtained for this unique 
system. Earlier interpretations were hampered by the lack of cer- 
tain vital information but the new X-shooter data have now al- 
lowed us a more complete view of this system. 

5. 7 . Geometry of the system 

Fig.Qsummarizes what we know about the redshifts of the com- 
ponents. The issue with the redshifts of the QSOs has been re- 
solved causing qA to be moved 740 km s"' towards higher red- 
shift. This means that qA is now offset only 100 km s ' from 
the Lyo- emission (at the position of qA) in the NOT spectrum, 
and 200 km s"' in the X-shooter spectrum. The difference is pre- 
sumably arising from the different way the velocity field of the 
emitter has been averaged over the slit in the two observations. 
The velocity off'set between the DLA and the Lyo- blob is much 
larger (640 + 70 km s '). All taken together the logical conclu- 
sion is that the blob is associated with qA and its host galaxy 
rather than with the DLA as it was previously thought. The end- 
to-end size of the blob along it's major axis is about 5 arcseconds 
corresponding to 42 kpc, and 28 kpc along its minor axis while 
the projected distance from qA to qB is 27.7 kpc. Based on the 
ionization state of the gas we find that the DLA cloud is placed 
along our line of sight towards qA at a distance larger than ~ 30 
kpc. 

5.2. Dynamical state of the system 

The systemic redshift difference between qA and qB is very 
small (60 km s"' and consistent with zero to within less than 
2cr) suggesting that their relative movement is in the plane of 
the sky. This would imply that their physical separation is the 
same as their projected separation. The Lya blob has a veloc- 
ity gradient along its major axis and through qA. This may be 
caused by infalling gas which is being ionized where it falls into 
the ionizing cones of q A as it was described in the models by 
IWeidinger et al.1 (l2004 l2005h but it could also be the signature 
of a rotating disk. From Fig. [3] and considering the resolution 
of the data along the slit corresponding to a smoothing length of 
13 kpc, we see that the current data are indeed consistent with 
the rotation curve of a large spiral galaxy in todays universe. A 
spectrum with better spatial resolution, and along the major axis 
of the blob, would be needed to discriminate between the two 
possibilities. 

The DLA absorber is a separate object which is moving to- 
wards the Lye blob with a relative velocity of 640 + 70 km s"' . 
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Fig. 10. Density profile of the two DM halos around qA and qB 
assuming that they have standard NFW profiles. 



5.3. Black Hole masses, DM Halo masses and Eddington 
ratios 

The mass of the DM halo of the QSO Q1205-30 at z = 3.04 
was determined from a model of the gas falling into its potential 
well and was found to be in the range log MtotlM^ - 12.3 - 12.8 
dWeidinger et alJ2004l) . As discussed in 35.2l it is not clear in our 
case if the dynamical signature of the gas is caused by infall or 
rotation so we shall use a different method to determine the DM 
halo mass. 

First we derive Black hole (BH) masses following 
IVestergaard & PetersonI (l2006l formulae 5 & 7) and use the re- 
lation between BH mass, HjS and Civ emission lines width 
and the continuum luminosity at 5100 A and 1450 A rest wave- 
length. We also estimated BH mass using the relation be- 
tween Mgii and continuum luminosity at rest-frame 3000 A 
(IVestergaard & Osmeill2009l formula 1). For the line widths we 
use the FWHM of the fitted gaussian for each emission line and 
we list those in TablelH We correc ted the FWHM of each Une for 
spectral resolution as described in lPeterson et al.l (l2004l) . For qB 
we do not have the H^S emission line width, therefore, we have 
no mass estimation from that line. 

To obtain the continuum luminosity at 5100, 3000 and 
1450 A we u s e the / /-band magnitudes of the two QSOs from 
Warren et al. (l200ll) and the / and B-band magnitudes from 



Fvnbo et al. (I2000h . We also checked for line contamination in 



these band and for this w e used the composite spectrum from 
IVanden Berk et al" I (l200lh . We subtracted emission line compo- 
nents from the spectrum and found the contamination to be 14, 
15 and 7% in the H, I and B bands respectively. After subtract- 
ing the contamination the resulting BH masses are listed in Table 
m It has recently been established that BH mass correlates well 
with the underlying total gravitational ma ss (presumably domi - 
nated by the dark matter) of the host halo jBandara et al.ll20()9l) . 
We use this correlation to finally obtain the corresponding halo 
mass of each QSO by using their mean BH mass and those are 
also listed in Table H] 

In Fig. [To] we plot the density pr ofile of the two halo s as- 
suming that they have NFW profiles jNavarro et alJI 19971) and 
that their projected distance (27.7 kpc) is also their real proper 
distance as argued in ^5.21 It is seen that qB is inside the halo 
of qA, but that it is still clearly defined as a separate sub-halo . 
We used the prescription from iNemmen" & BrothertonI (1201 Ol) 



and iGavignaud et al.l (l2008h to determine bolometric luminosi- 
ties and eddington ratios (provided in Table |4ll respectively. 



6. Discussion and Conclusions 

In this paper we have presented spectra of the two QSOs 
Q 015 14-048 A and Q0151-H048B as well as of a Lya blob at 
the same position. We obtained spectra of the blob at two dif- 
ferent position angles. QSO redshifts determined from Hy6 and 
Ha change the previously reported redshifts by 740 km s ' in 
ad dition to th e much smaller corrections previously suggested 
by iTvtler & Fa n (1992). Our first conclusion is that the correc- 
tions to QSO redshifts based on rest-frame UV lines can be much 
larger than previously thought and that other QSO absorption 
systems of this kind (i.e. with large inverted redshift offsets) pre- 
sumably will turn out to find the same resolution. 

The new redshift of qA means that the Lye blob is now 
more likely belonging to this QSO. We determine a Black Hole 
mass of qA of 10^ "'^ Mq and a corresponding DM halo mass of 
10'^^"^ Mq. The Lyor blob is sitting in the potential well of this 
halo. The blob has a size of 42 by 28 kpc, it has a clear signature 
of a strong velocity gradient along its major axis which is con- 
sistent with the rotation curve of a large disk galaxy in todays 
universe but may alternatively be caused by in- or out-flow. In 
the X-shooter spectrum along the minor axis which has higher 
spatial resolution there is no velocity gradient of the blob, but 
there is a hint that the Lya emission may be double peaked. 

Interpolating in the merger tree simulations of 
[Parkinson et alj ( l2008h we find that the halo of qA in the 
mean will end up in a halo with a mass of 6 x lO''* M© at z = 0. 
We report bolometric accretion rates of both quasars and find 
that they accrete at 60% and 30% of their Eddington limits. This 
places both qA and qB consis tently on the relation shown in 
fig. 4 of IGavignaud et al.1 (l2008h . 

From the aspect of the two quasars they therefore appear to 
be very typical for quasars at that redshift, and their environ- 
ment is apparently a proto-typical pre-cursor of what will be- 
come something very close to our local super-cluster at z = 0. 
What makes this system special is the accidental alignment of 
many different types of objects so close together The physi- 
cal interactions between the dual quasar, the Lyo- blob, the host 
galaxies and the PDLA absorber makes this a unique laboratory 
to study what seems to be what our local super-cluster looked 
like aiz-2. 

The metallicity of the DLA inferred from the X-shooter 
medium resolution spectrum is [Si/H] - -1.93 + 0.04, in 
the lowest 15% of DLAs at this redshift. This stands in con- 
trast to the recent result by Ellison et al. (2010) who reported 
that PDLAs in general have higher metallicity than intervening 
DLAs. 

Based on the ionization state of the gas we model the phys- 
ical state of the DLA and find a best fit at a distance of 30-40 
kpc from qA assuming that contribution of local star-formation 
can be igno red. Under the same assumption, and using a differ- 
ent method, lElUson etal J (120101) found that the distance must be 
larger than 160.5 kpc. In case there is local star formation close 
to the DLA our model would only provide a lower lirn it and 
would be consistent with the result of lEUison et alj ( l2010l) . 

The D LA falls nicely on t he relation of velocity-width vs 
metallicitv lLedoux et aP (l2006l) . Using the [X/H] vs the TJ-band 
magnitude relation derived in that paper we predict that the DLA 
should have an R-band brightness of about rriR - 30 and a Lye 
flux of around 2 x 10"'^ erg s"' cm"^. This object is therefore 



Tayyaba Zafar et al.: Q0151+048 system at z = 1.93 



11 



not observable in emission in our data, consistent with the inter- 
pretation of the Lya blob belonging to the host galaxy and not to 

the PL A. 

iHeckman et alJ (Il99lh reported the detection of "fuzz" 
around radio loud quasars and argued that the fuzz had to 
be clumpy on s mall s cales with a very small filling factor. 
IWeidinger et al.l (l2004l) reported similar fuzz around a radio 
quiet quasar and showed that its velocity profile and morphology 
implied that it was in a state of infall. That quasar also exhibit a 
strong absorber close to its emission redshift, a system which in 
current terminology may be named a proximate sub-DLA. The 
du mpiness on sma l l scale s and the small filling factor predicted 
by iHeckman et al. J 199lh matches the few known PDLA sys- 
tems well (qA, Leibundgut & RobertsonI 19991 Weidinger et al] 



12004 iHennawi et al.ll20()9ir and like IHennawi et al.l ( l2009l " we 
consider it likely that PDLA Lya blobs in some ways are related 
to the fuzz around radio loud quasars. However, there are signif- 
icant differences in terms of energy output, in the measured Lya 
Une widths (which is much narrower for qA), in the mechanical 
interaction with the radio jet and in the fact that metal emission 
lines are detected in the fuzz around rad io loud object. The case 
of qA is therefore more similar to that o f lWeidin^er etaH (l2004 
and the negative velocity of the DLA is an indication that it is in 
a similar state of infall. 

We detect emission in the trough of the DLA which we in- 
terpret as the signature of the quasar host galaxy. It is interesting 
to note that the hosts were not detected in any of the three above 
mentioned cases which are all at z = 3. This is consistent with 
recent reports that there is a rapid transi tion of Lyg selected ob- 
jects o bserved at a redshift of about 2.5 (iNilsson & M0lledl2OO9l 
I20TI . 

We briefly report the detection of four intervening C iv sys- 
tems seen in both sightlines of which one system has Si iv ab- 
sorption seen only in a single sighline. Also we see rich Lyman 
forests with near identical absorption signatures. This QSO pair 
is therefore well suited for a multi-sightline study of both Lyman 
forest absorbers and of intervening metal absorbers. 
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